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Abstract 

We present a calculation of the complete O (a) electroweak radiative corrections 
to the M0ller scattering process e~e~ — > e~e~ , including hard bremsstrahlung contri- 
butions. We study the effects of these corrections on both the total cross section and 
polarization asymmetry measured in low energy fixed target experiments. Numerical 
results are presented for the experimental cuts relevant for E-158, a fixed target e~e~ 
experiment being performed at SLAC; the effect of hard bremsstrahlung is to shift the 
measured polarization asymmetry by ~ +4%. We briefly discuss the remaining theo- 
retical uncertainty in the prediction for the low energy M0ller scattering polarization 
asymmetry. 
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1 Introduction 



The search for physics beyond the Standard Model (SM) can be pursued in two distinct 
ways: by observing the direct production of particles associated with this new physics, or 
by detecting the indirect influence of these new states on precision measurements. Since the 
discovery of the top quark in the early 1990s, no new elementary particle has been found. 
However, many indirect constraints on the SM and its possible extensions have been obtained 
by experiment. For example, high precision measurements on the Z-pole at LEP, combined 
with the values of My/ and m t measured at the Tevatron, have constrained the mass of the 
SM Higgs boson to satisfy the bound ran < 196 GeV at 95% CL by studying its influence on 
the radiative corrections to Z — > // decays pQ. Similarly, the parameter space available to 
super symmetric theories is becoming restricted by the absence of CP-violating electric dipole 
moments which are required in supersymmetric models to explain the baryon asymmetry of 
the universe via electroweak baryogenesis E] ■ 

There are currently several quantities whose measured values do not match their SM 
predictions. The value of the muon g - 2 |j, sin 2 (9y/) as derived from the low energy 
neutrino- nucleon scattering experiment NuTeV [Hj, and the magnitude of atomic parity vi- 
olation (APV) in the cesium atom [Hj have all been recently reported to deviate from their 
SM prediction by > 2.5 a. The interpretation of these discrepancies, however, is rendered 
difficult by our current inability to reliably calculate the required hadronic or atomic physics 
that affects these quantities. The deviations discussed above have been argued to either be 
less significant or to vanish altogether when the relevant contributions are more carefully 
analyzed [7]. 

These results increase the importance of E-158, a fixed target e~e~ experiment cur- 
rently being performed at SLAC [8J. This experiment will determine sin 2 (9y/) at a mo- 
mentum transfer Q 2 « 0.02 GeV 2 , comparable to the value relevant for neutrino-nucleon 
scattering, by measuring the polarization asymmetry A LR = (o L — a R ) / (<j l + a R ) in the 
M0ller scattering process e~e~ — *> e~e~; this reaction is less sensitive to the hadronic 
physics that complicates the previous measurement of sin 2 (0w) by NuTeV, and its inter- 
pretation should therefore not be plagued by uncontrollable theoretical uncertainties. The 
error on the E-158 measurement is expected to reach 5Alr/Alr « ±8%, corresponding 
to 5 sin 2 {6 w) = ±0.0008, making it the most accurate determination of sin 2 (8w) at low 
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momentum transfer 

The study of the electroweak (EW) radiative corrections to M0ller scattering was 
initiated by Czarnecki and Marciano jTO]. They found that the one-loop corrections reduce 
the tree-level prediction for Alr by « 40% at low Q 2 in a MS renormalization scheme; the 
large size of this effect can be traced to the fact that while the tree level prediction is pro- 
portional to the numerically small electron vector coupling, g v = —1/4 + sin 2 (9w) ~ — 10 -2 , 
the one-loop result contains quark contributions to 7 — Z mixing not similarly suppressed. 
Although the light quark components cannot be computed in perturbation theory, they can 
be related to e + e~ — > hadrons scattering data; the analysis of JU] concludes that these effects 
can be determined with adequate precision. The virtual electroweak corrections to M0ller 
scattering for arbitrary energies were calculated in ^T], while the pure QED component of 
e~e~ — > e~e~ was computed in ^2] without experimental cuts. Neither of these results is 
sufficient for comparison with the E-158 measurement. 

In this paper we compute the complete O (a) EW corrections to M0ller scattering, 
including hard bremsstrahlung effects, and impose the experimental cuts relevant for the 
E-158 measurement at SLAC. Our recalculation of the virtual EW corrections is valid for 
arbitrary center-of-mass energies and serves as a check of the results in [11] . with which we 
agree. We find that the majority of the QED corrections to Am cancel when both virtual 
and hard photon corrections are added, as anticipated in [TUl HI] ; they contribute only a 
small ~ +4% shift to A LR . This is due to the experimental setup of E-158, which treats 
photons inclusively. 

The paper is organized as follows. In Section 2 we present our notation and briefly 
discuss the tree-level M0ller scattering process. We describe the calculation of the one-loop 
electroweak corrections in Section 3, and explain the treatment of the hard bremsstrahlung 
effects in Section 4. We present the appropriate experimental cuts and numerical results in 
Section 5, and conclude in Section 6. 

2 Preliminaries and leading order results 

In this section we present the Born-level predictions for the unpolarized cross section and po- 
larization asymmetry for the process e~(pi)e~(p 2 ) ~^ e ~(<li) e ~(Q2)- The notation introduced 
here will be used throughout the paper. 
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We first introduce the following Mandelstam invariants: 



s = (pi+p 2 ) 2 , s' = (q 1 + q 2 ) 2 , 
t = (pi ~ qif , t = (pa - q 2 ) 2 , 

u = (pi - g 2 ) 2 , u' = (P2 - <?i) 2 • (1) 

In the 2 — >• 2 scattering process described by the Born-level prediction and the virtual 
corrections, the primed invariants are equivalent to their unprimed counterparts: s = s , 
t = t , and u = u . These relations will not hold when we consider the corrections arising 
from photon emission. We next introduce the electron vector and axial couplings to the 
Z boson, g v = — 1/4 + Syy and g a = 1/4, where = sin 2 (0w) is defined in the on-shell 
renormalization scheme: = 1 — M^/Mf. It is convenient to define the following energy- 
dependent effective couplings: 



CLL[X} *\x Jt sicl U-.V- 



1 {9v ~ 9a 



,2 
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crr (x) = e 2 



1 (9v + 9a) 



clr(x) = e 2 {-+ 2 £> ) , (2) 

[x s 2 w c 2 w (x - M|) J 

where — e is the charge of the electron. Using these expressions, we can write the squared 
tree-level matrix elements as 

\M° LL \ 2 = A{c LL {t) + c LL {u)) 2 s 2 , 
\M RR \ 2 = 4(c RR (t) + c RR (u)) 2 s 2 , 

\M° LR \ 2 = \M RL \ 2 = A(c LR (t)) 2 u 2 + A(c LR (u)) 2 t 2 . (3) 

The subscripts encode the polarizations of the initial electrons, with L denoting left-handed 
states and R representing right-handed states. All terms of Oiw? e jx\ where m e is the 
electron mass and x denotes one of the invariants introduced in Eq. ^ have been neglected. 
This is an excellent approximation in the relativistic limit, and is appropriate for the E-158 
experimental setup. We next define polarized differential cross sections, 
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(4) 



and let denote the integration of these expressions over the region under consideration. 
We remind the reader that the cross section is invariant under boosts along the beam axis if 
the boundaries of integration are chosen consistently in each frame. The unpolarized cross 
section and polarization asymmetry can now be written as 

1 r i 

T {&LL + &LR + &RL + CTRRf > 

&LL + &LR — &RL — &RR , rN 

• ( 5 ) 

&LL + 0~LR + Cr RL + a RR 

Before discussing the electroweak corrections to these quantities, we comment briefly 
on our on-shell choice of s^. Numerically, ~ 0.2216 in the on-shell scheme, while the 
MS value is (Mz)jjs ~ 0.2312. Recalling that the tree-level A LR is proportional to 
g v = — 1/4 + sin 2 (9w), we see that the choice of the on-shell renormalization scheme tends 
to increase the tree- level asymmetry; the value of A LR is also affected by our choice of M w 
as an input, rather than G M as in JU] (this choice is made to facilitate comparison with the 
complete EW virtual corrections computed in jTT]; the additional uncertainty induced by 
the larger experimental error in My/ is small compared to the uncertainty arising from the 
hadronic contributions to 7 — Z mixing). We will find that the electroweak corrections to 
the Born-level prediction are larger in the on-shell scheme than in the MS scheme chosen 
in pH] (this enhancement of on-shell radiative corrections over those of the MS scheme is 
familiar from other studies of EW physics JSj)- This indicates that the prescription of JU] 
is probably a more appropriate parameterization of the tree-level prediction; however, since 
we will primarily study the effects of hard bremsstrahlung on A^ R , this fact will not concern 
us here. 

3 Virtual corrections 

In this section we study the one-loop electroweak corrections to the M0ller scattering process. 
We reduce all one-loop integrals to the scalar Passarino-Veltman basis using FORM |15j : 
the numerical evaluation of the resulting integrals is then performed using LoopTools |16j . 
The matrix elements obtained are valid for yfs ^> m e . Since the resulting expressions are 
rather lengthy, and have been previously discussed in the literature jllj . we will not present 
them explicitly in this manuscript; they are available upon request from the author. We 
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Figure 1: Representative diagrams contributing to the one-loop QED corrections to M0ller 
scattering. 

use dimensional regularization to regulate ultraviolet divergences, and an anticommuting 75; 
this is appropriate since no ABJ anomalies are present |17j . 



3.1 QED corrections 

We follow fl] and define the QED corrections to consist of those diagrams where one photon 
leg has been attached to a tree-level diagram. A representative sample of the relevant graphs 
is presented in Fig. ^ We work in the complete on-shell renormalization scheme of ^H] , in 
which self-energy corrections to external fermions do not explicitly appear. 

These diagrams contain both infrared (IR) and ultraviolet (UV) divergences. We regu- 
late the IR divergences with a finite photon mass \i\ these cancel against identical divergences 
appearing in the process e~er — > e~e~7, where the photon has an energy k < AE <C y/s. 
In this soft photon limit, the integral over the photon phase space factorizes, and its effect 
is to multiply the tree- level squared matrix elements by the following correction factor A s : 

\M s t] \ 2 = A S \M%\\ (6) 

where A s has been derived in ^T] and is given by 

The appearance of the cutoff energy AE indicates that this expression is frame-dependent. 
The E-158 experiment measures the fully inclusive cross section e~e~ — *> e~e~ + (n)j; the 
dependence on AE will therefore cancel from our results when hard bremsstrahlung effects 
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are included, and we need not consider this issue. The UV divergences can be cancelled 
by the following wavefunction renormalization of the electron field: ip e — ► (1 + 1/2 8Z^) ip e - 
This renormalization introduces counterterms for both the jff and Zff vertices. In the 
on-shell scheme, SZ^ can be written in terms of Passarino-Veltman (PV) functions as 

5Z^ = -^{B (0,m e ,0)-l-Am 2 e B' Q {mlm e ,fj}j} , (8) 

where a is the fine structure constant, B is the two-point PV function in the notation of |18j . 
and B Q is the derivative of this function with respect to its first argument. B is infrared 
divergent, and we have therefore retained its dependence on the photon mass \x. 

Denoting the one- loop QED matrix elements by M-fj, we can write the QED correc- 
tions to l-Mjjl 2 as 

5f ] =2Re{Mf J (MlY} + \Mt,\ 2 . (9) 



3.2 Weak corrections 



The one-loop weak corrections to M0ller scattering consist of the 7 and Z self-energy graphs, 
7 — Z mixing, the remaining vertex and box diagrams, and the appropriate self-energy and 
vertex counterterms; a representative subset of the countributing graphs is presented in 
Fig. |21 These contributions contain UV divergences, which must be removed via the full 
SM renormalization program; we follow the on-shell prescription described in JH]- The self- 
energy and counterterm graphs can be combined and interfered with the Born-level diagrams 
to yield 



M° LL (M S L + L C )' 



M° LR (Ml+ R c 



4e 2 ({c LL (t) + c LL (u)} {V LL (t) + V LL («)}) s 2 , 

4e 2 ({c RR (t) + c RR (u)} {V RR (t) + V RR (u)}) s 2 , 

M° RL (M S R + L C Y = 4e 2 (c LR (t) V LR (t) u 2 + c LR (it) V LR (u) t 2 ) , (10) 



where the (x) are given by 
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Figure 2: Representative diagrams contributing to the one-loop weak corrections to M0ller 
scattering. 
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(11) 



We have introduced the abbreviations ^ = (g v — g a ) /syycw and g# = (g„ + g> ) / sy/Cw for 
the left and right-handed electron couplings. The forms of the self-energy insertions S AA (x) , 
Yj AZ (x), and Y, zz (x), as well as expressions for the counterterms 5Z e , SZl, SZ r , 5sw, and 
5 Ml in terms of these functions, can be found in ^H]; note that the pure QED component 
of Eq. ^2 has been considered separately in the previous subsection. 

The momentum transfers relevant for the E-158 measurement are \t\, \u\ ~ 0.02 GeV 2 . 
At this energy scale, the light quark contributions to the 77 and 7Z vacuum polarization 
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functions in Eq. [TT] cannot be computed perturbatively; the appropriate degrees of freedom 
in this regime are the low-lying hadronic states, and the required terms must be obtained by 
comparison to experimental data. We discuss here the determination of these contributions. 
We first consider the expression 2 5Z e — T, AA (x) jx. Evaluating the hadronic contributions 
to this quantity in the free quark approximation, and neglecting the quark masses wherever 
possible, we find 

»*-^-^-{ ta S)-5} + -- (12) 

where the sum is over the five light quarks and the ellipsis denotes the terms arising from the 
remaining states. This quantity can be recognized as A had (x) = c%ad (0) — a^ad which 
can be related to e + e~ — > hadrons scattering data via a dispersion relation [TH] . We use the 
parameterization of Ahad (%) given in [20 j , which is a good approximation throughout the 
entire range of x. A similar calculation reveals that 

2 5Z e A had (Mf) + E AA ( Mf ) I M\ + 



t - Ml t-Ml 

where only the light quark components have been explicitly shown; Ah a d iM"z) can again be 
obtained from [H|, while Z AA (Mf) can be calculated perturbatively. 

Finally, we must discuss the contribution of the 7— Z mixing term, — £ (x) jx {x — Ml] 
For \x\ much larger than the light quark masses, this quantity can be evaluated perturba- 
tively; as mentioned earlier, this condition does not hold in the E-158 experimental setup. 
We begin by rewriting this piece using 



E AZ (x) 1 



x e 



+ I = !&M + I, (i4) 

e x e 



where e = 4 — D is the usual regulator appearing in dimensional regularization, and U.]^ s 



is the 7Z vacuum polarization function defined in the MS renormalization scheme. The 
remaining 1/e pole cancels when the counterterms of Eq. are combined with the vertex 

diagram contributions. We next set x = in II jj^; an estimate of the error induced by 
evaluating this quantity at x = rather than at the value \x\ ~ 0.02 GeV 2 relevant for the 
E-158 experiment was performed in JUj by calculating the pion contributions to IP Z (x) in 
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a scalar QED framework. This effect was found to be negligible, and we ignore it in the 
remainder of our analysis. Evaluating IP Z in the free quark approximation, we obtain 



N c Q q T 3 -2Q q s 



2 

q°W 



3 



Ml 



(15) 



A dispersion relation analysis of this quantity was performed in [TTH I21j; we use this result 
and replace 



\j2 N cQi [T 3 -2Q q s 

6 Q 



Ml 



-6.88 ±0.50 



(16) 



It is claimed in that an improved treatment of the e + e _ — > hadrons scattering data 
would reduce the quoted error; an analysis of this type motivated by the needs of a high 
energy e~e~ program seems to be underway 



Denoting the one- loop weak matrix elements by , we can write the weak correc- 



tions to | .M^ | 2 as 

< = 2Re {.<(<)*} . 
The complete one-loop expressions for the matrix elements become 



(17) 



\M 



°.\ 2 + s9. 



Our results are in a form valid for all ^/s ^> m e ; we can compare them with the detailed 
results for the unpolarized cross sections and polarization asymmetries given in Tables 1 
and 2 of jllj . When we use the same parameters found there, the older parameterization of 
Ahad (x) given in , and the light quark masses given in ^1 > we find complete numerical 
agreement to the given accuracy. 



4 Hard bremsstrahlung corrections 

In this section we describe corrections to M0ller scattering arising from the process e~(pi)e~(p2) 
e ~{ ( li) e ~{ < l2)l (k), in which the emitted photon has an energy k > AE. We first discuss the 
calculation of the relevant matrix elements, and afterwards present our parameterization of 
the phase space. 
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In our calculation of the hard bremsstrahlung corrections to M0ller scattering, we 
again work in the relativistic limit y/s ^> m e ; however, as is well known, m e cannot be 
neglected in regions of the phase space where the emitted photon becomes collinear with one 
of the electrons. Similarly, m e must be retained when one of the final-state electrons travels 
parallel to the beam axis, and the t-channel momentum transfer becomes of O (m 2 e ). We first 
perform the calculation neglecting m e everywhere, and then discuss the required correction 
factors. We introduce the following abbreviations for the electron propagator denominators 
connected to an emitted photon: 

-1 -1 
di = , d- 



2pi ■ k 2p2 ■ k 

' <k = —— . (19) 



2q 1 ■ k 2q 2 ■ k 

The expressions for the squared matrix elements |A^' a | 2 computed with m e = everywhere 
can be found in Eq. ESI in the Appendix. Recalling the definition of (x) presented in 
Eq. El we observe that the matrix elements contain peaks of the form s/m 2 e when one of 
the following eight quantities becomes of O (ml): l/di,t,u,t' ,u , where i = 1,2,3,4. The 
first four peaking regions are associated with the emitted photon becoming collinear with 
any of the four electrons, while the remaining four peaks arise when one of the final-state 
electrons becomes collinear with an initial-state electron. Note that all terms containing 
d 2 ,l/t 2 ,l/u 2 ,l/t' 2 , and 1/u 2 have cancelled, leaving no contribution which behaves like 
s 2 jm\. However, in the peaking regions, terms of the form m 2 e d 2 , m 2 e /t 2 , etc. are of O (s/m 2 e ), 
and must be included; we now discuss the calculation of these pieces. 

We first consider the mass effects associated with the unpolarized final-state particles. 
We focus on e~ (gi); the discussion proceeds identically for e~ (q^)- A typical final-state 
bremsstrahlung diagram is shown in Fig. 01 the self-interference of this diagram, or its 
interference with another diagram where 7 (k) is emitted from e~ (qi), leads to the following 
term in the squared matrix element: 

{;/!+ /^ + m e } ^{k)u{q 1 )u{q x ) f (k) {^+ $ + m e } (d 3 ) 2 
= {2 Ql -e(k)+ It f(k)}u(q 1 )u(q 1 ){2q 1 -e*(k)+ ft f (k)} (d 3 ) 2 , (20) 

where e(k) is the photon polarization vector, u (q%) is the Dirac spinor of e~ (qi), and the 
equality results from the use of the Dirac equation (^1 + m e ) u (qi) = 0. We now sum over 
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7(A) 



e (pi) / e (qi) 



e (p 2 ) / e (q 2 ) 



Figure 3: A typical bremsstrahlung diagram; the self-interference of this particular graph 
neccessitates the inclusion of final-state mass effects. 

photon polarization vectors and electron spin states; since the electron is unpolarized, we 
replace u(q{)u (gi) — >q{ + m e . After standard Dirac algebraic manipulations, keeping only 
l/c?3 and m 2 / (d^) 2 terms, we arrive at the expression 

2 &d 3 -4m*(#+ ^)(d 3 ) 2 . (21) 

We have obtained the required l/ds and m 2 / (d^) 2 terms; the remainder of the numerator 
algebra can now be performed in the massless limit. We note that d% must be evaluated 
with its full m e dependence in the peaking region; otherwise, d% — > oo, and an unphysical 
singularity is encountered. 

The calculation of initial-state mass effects is similar to the process outlined above. 
We find that the squared matrix element contains the term 

{2 Pl ■ e (k) - # /(k)} u ( Pl ) u (pO {2 Pl ■ e* (k) - # f (k)} (d,) 2 (22) 

when self-interference of diagrams with bremsstrahlung off e~ (pi) is considered; a similar 
expression is obtained when studying e~ (p 2 ). Since the initial states can be polarized, we 
must now replace 

u (pi) u (p^ -> ^ (fa + m e ) (1 + 75 f) , (23) 



where 

r 



P|, ^-mA , (24) 
m e Pi • ? I 



P\\ is the degree of longitudinal polarization, and q is an auxiliary vector satisfying q 2 = 0, 
Pi • g t^O; for mass effects associated with e~ (pi), it is convenient to choose q = p 2 . For the 
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calculation of mass effects arising from e~ (P2), we choose q = pi- After performing standard 
Dirac algebraic manipulations, again keeping only the leading terms which contribute to the 
squared matrix element, and setting Pu = ±1, we find that Eq. l22l becomes 

A 2 

u± {_2 fa -4m 2 (^- 16) (d 1 ) 2 } T ^(u+- w_) { Pl ■ p 2 ft - p 2 ■ k #)} (d,) 2 , 

(25) 

where u± = (1 ± 75) /2 are the spin projection operators. We find a similar expression 
when calculating the mass effects arising from e~ (P2). The explicit expressions for the mass 
correction terms |.M^ ! | 2 can be found in Eq. EH] in the Appendix. 

Finally, we must discuss the mass effects which are relevant when t, u, t , or u becomes 
of 0(m 2 ). Such terms arise only from photon exchange diagrams, and cancel from the 
numerator of Alr. We consequently need not worry about any complications arising from 
polarized initial beams. Eq. already contains terms of the form m^c^ (x) d\ — ► mldl/x 2 , 
which become important when either x or 1/dk is of O (m 2 ). We need therefore only derive 
the terms of the form m 2 didj/x 2 . It is possible to have "overlapping peaks", where both x 
and 1/di become of O {m 2 e ) at the same point in phase space. Such events are characterized 
by a photon traveling in one direction down the beam axis while an electron goes in the 
opposite direction; momentum conservation then implies that the remaining electron has 
a very low energy. The experimental cuts relevant for the E-158 measurement require the 
detection of at least one electron with a reasonably high energy at a large center-of-mass 
frame angle, thereby eliminating all overlapping peaks. The appropriate mass terms to add 

to the total unpolarized cross section, |.M^ pol | 2 , can be found in Eq. 1571 in the Appendix; 
the same contribution enters the denominator of A^r. 

We now have all of the pieces needed for the hard bremsstrahlung cross section. The 
matrix elements appearing in the numerator of the polarization asymmetry are 

IM'l 2 = \Mlf\ 2 + \Mlf\ 2 ; (26) 

the matrix elements needed for the unpolarized cross section and the asymmetry denominator 
become 

\M h unpol \ 2 = J2{\ M if\ 2 + \Mlf\ 2 } + |AC pol | 2 . (27) 

y 

We must now present our parameterization of the e~e~ — > e~e~7 phase space. We 
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follow [21] and write the hard bremsstrahlung cross section for a given helicity configuration 

as 



(28) 



where g = gi + g 2 and the two sets of integrations are separately Lorentz invariant. We 
choose to evaluate the left set in the center-of-mass (CM) frame, and the right set in the rest 

frame of the qi + g 2 system. Doing so, using s = (gi + g 2 ) 2 = [p\ + p 2 — k) 2 = s — 2E k ^, 
and performing an azimuthal integration, we arrive at 



4 = - - , 4 / dE k dc k dQ* E k \M h i3 \ 2 . (29) 



1 

64s (2tt) 4 

Q± 2 * s the rest-frame solid angle of the final-state electron e~ (gi); all rest-frame quantities 
will be denoted by a superscript R, while all lab (fixed target) frame quantities will have 
a superscript L. Expressions without superscripts, such as the photon energy E k or the 
cosine of its production angle c k = cos(9 k ), denote CM frame quantities. In the absence of 
experimental cuts, the limits of integration are 

s — Ami 

AE < E k < 6 , -1 < c k < 1 , 

-1 < c f = cos (Of) < 1 , < 0f < 2tt . (30) 

The relevant cuts involve restrictions on the lab frame electron energies and angles; we 
implement them numerically. It is a simple task to express the Mandelstam invariants of 
Eq. n m terms of these integration variables, and to derive the Lorentz transformations 
connecting the lab, CM, and rest frames. 

The complete expressions for the unpolarized cross section and polarization asymme- 
try become 



On 
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A a LL + O h LL + °LR + gig ~ Q RL ~ Orl ~ O KR ~ Vrr 

A LR - , (Olj 

4cr„ 
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where o^- is the 2^2 cross section defined in Eq.|H and the one-loop matrix elements |-Mr,-| 2 
are defined in Eq. We can now study the effects of both virtual and hard bremsstrahlung 
corrections on measurements of the total cross section and polarization asymmetry in low 
energy M0ller scattering. 

5 Numerical results 

In this section we present numerical results for the experimental setup relevant for the E-158 
measurement at SLAC. For a detailed description of the E-158 experiment, we refer the 
reader to [HI E] ; we give below a description of the cuts relevant for our purposes. We remind 
the reader that we have split the radiative corrections into QED and weak components; we 
will study the effects of these two types of corrections separately. 

We demand that at least one electron be detected, with its lab frame energy and 
production angle satisfying the following constraints: 

E L > 10 GeV , 4.65 mrad < 9 L < 8.00 mrad . (32) 

The exact form of the angular cut is chosen so that typical scattered electrons have a CM 
frame angle in the range —0.5 < c\ < 0. The structure of the E-158 detector is such that 
events where both scattered electrons enter the acceptance defined by Eq. E21 are counted 
twice; with this restriction on Q\ , no 2 — > 2 scattering process leads to such an overcounting. 
However, events with hard photon emission can create such an effect; we will study their 
contribution to both the total cross section and A LR . We choose a beam energy E b = 48 GeV, 
and the following gauge and Higgs boson masses: M w = 80.451 GeV, M z = 91.1875 GeV, 
and itih = 120 GeV. We use the fermion masses and fine structure constant given in (22] • 
With this choice of parameters and cuts, we find the following results for the Born-level and 
weak-corrected unpolarized cross section and polarization asymmetry: 

al = 15.348 /ib , A B LR = 3.646 x 10~ 7 , 
a B+w = 15.589 nb , Af+ W = 1.756 x 10~ 7 . (33) 

The effect of weak radiative corrections is to decrease the polarization asymmetry by ~ 52%. 
A 48% decrease from weak corrections was found in ^1] using the on-shell scheme; this result 
is for the specific phase space point c\ = 0, different values for Mw, Mz, and mn, and an 
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older parameterization of Aahad- The two results agree when these differences are taken into 
account. We will compare our results with the MS calculation of [10J after discussing the 
QED corrections. We note that the dependence of both the total cross section and Alr on 
ran is weak; choosing mu = 200 GeV leads to an « 3% increase in Alr and a neglgible 
change in the cross section. 

We present in Table [0 the QED-corrected cross section and asymmetry, and show 
the percent shifts from the weak-corrected values (J^ +w and Afj^ given in Eq. EH1 Both 
virtual and hard bremsstrahlung corrections are included. In our bremsstrahlung calculation 
we introduced an energy cutoff AE, below which the emitted photon phase space integrals 
were performed analytically in a soft-photon approximation, and above which the exact 
bremsstrahlung matrix elements were treated numerically. Since we integrate inclusively 
over the photon phase space, the dependence on AE should cancel; we include in Table [Hour 
results for four choices of this parameter to verify that this indeed occurs. These results were 



AE (yTs) 


a u (/lb) 




A LR (10- 7 ) 


S Q A LR 


10~ 6 


15.455 


-0.86% 


1.817 


+3.5% 


10~ 5 


15.456 


-0.85% 


1.816 


+3.4% 


10" 4 


15.457 


-0.85% 


1.817 


+3.5% 


10~ 3 


15.457 


-0.85% 


1.817 


+3.5% 



Table 1: Total unpolarized cross section, a u , and polarization asymmetry, Alr, for different 
values of the photon energy cutoff AE. Also included are the shifts in these quantities from 
(Tu +W an d Afft W , defined in Eq. induced by QED corrections. AE is given in units of 
y/s, and a u is given in millibarns. 

obtained using the Monte Carlo integration program VEGAS [26J, with the number of calls to 
the integrand N caU = 5 x 10 6 . The statistical error should therefore equal 1 / \J N caU « 0.05%; 
we have included enough digits in the values of Table ^ to test this expectation. We see 
that the variation of a u and Alr with AE is consistent with this statistical fluctuation. 
The shift in the total cross section induced by QED corrections is small and negative. The 
infrared-safe combination of virtual corrections and soft photon emission reduces the cross 
section, while the hard photon cross section partially compensates. These results are similar 
to those of [2Zj for Bhabha scattering, where total QED corrections to the cross section of 
~ —5% were found for —0.5 < c% < 0, with an acollinearity cut of 15° on the angle between 
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the final state electron and positron and an energy cut on the detected electrons. Our lack 
of an acollinearity cut allows kinematical configurations in which an undetected electron 
travels down the beam axis, and either u, t, u , or t becomes of O (ml), adding an additional 
positive contribution to the QED correction. For illustrative purposes we present the total 
cross section for three different acollinearity angles in Table 121 The decrease in the cross 
section is reduced as the acollinearity angle is increased, and a larger amount of the photon 
phase space is included. 



@acoll 


a u (lib) 


S Q a u 


5° 


13.865 


-11% 


10° 


14.407 


-7.6% 


15° 


14.688 


-5.8% 



Table 2: Total unpolarized cross section, a U} for different choices of the acollinearity angle 
Oacoii- Included are the percent shifts from the weak-corrected value in Eq. EHl 

We observe from Table Q that QED corrections increase the prediction for A^r by 
~ 3.5%, leading to a final prediction of A^r = 1.82 x 10~ 7 in the on-shell scheme. The 
one-loop MS result of |10] is A LR = (1.80 ± 0.09 ± 0.04) x 10~ 7 , where the quoted errors 

arise primarily from the hadronic uncertainties in W z \ this value does not include either 
bremsstrahlung effects or experimental cuts, and is evaluated at the specific phase space 
point cos (8) = where the tree-level prediction for Alr is maximized. The very close 
agreement between these results is surprising and to an extent accidental; however, we can 
identify the following two physics effects which contribute to the agreement: 1) the significant 
difference between the tree- level prediction for Alr in the on-shell and MS renormalization 
schemes caused by differing choices of sin 2 (9^) is ameliorated once virtual corrections are 
included, and 2) the experimental cuts relevant for the E-158 measurement render hard 
bremsstrahlung corrections a small effect. Also, the ~ +3.5% shift induced by hard photon 
emission accidentally brings the weak-corrected value of Alr in Eq. 133] into closer agreement 
with the result of 1 0J . QED corrections to the photon exchange component of M0ller 
scattering of ~ 100% are reported in [T2j. We believe that the large size of these results 
arises from the lack of experimental cuts in this calculation; without such restrictions, large 
contributions from the peaking regions discussed in the previous section greatly enhance the 
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cross section. The effects observed in a realistic experimental setup are significantly smaller, 
as our calculation demonstrates. 

We now briefly review the error on the theoretical prediction for Alr- The uncertainty 
in the one-loop result arising from hadronic contributions to the 7Z vacuum polarization 
function has been mentioned above and discussed in the previous section; although large, 
there is hope that it will be reduced in the future j22j. The remaining errors can be divided 
into those arising from either higher-order QED or weak corrections. The small magnitude 
of the QED corrections found in this paper, « 3.5%, indicate that the inclusion of O (a 2 ) 
effects are unnecessary; the experimental cuts protect the E-158 measurement from large 
QED effects. The one-loop weak corrections are large: ~ 40% in the MS scheme and « 50% 
in the on-shell scheme. However, the size of these effects arises from quark contributions 
to Alr not suppressed by 1 — 4 sin 2 (6w), a qualitatively new feature which first appears at 
0(a). No similar effect appears at the next order in perturbation theory, and we expect 
that the higher order corrections are suppressed by a factor of a/n ~ 0.1% relative to 
the one-loop result^. We therefore conclude that the most significant source of theoretical 
uncertainty on the measurement of Alr arises from the hadronic contributions to the 7Z 
vacuum polarization. 

As mentioned above, the E-158 measurement double-counts events where a hard 
photon knocks both electrons into the detector acceptance; this modifies the prediction for 
both a u and Alr. We find that this leads to the following effective values of the cross section 
and asymmetry, as well as the following shifts from the weak-corrected results of Eq. 1221 

a e J f = 15.613^b , 5f ff a u = +1.5% , 
A e H = 1.829 x 10- 7 , 6? ff A LR = +4.2% . (34) 

We have chosen AE = 10~ A y/s to obtain these values, and have checked that the variation 
with this parameter is consistent with the statistical error of the integration. The overcount- 
ing of events leads to an effective increase of the cross section, and increases the shift of 
Alr to +4.2%; subtracting from this the ~ 3.5% shift of Table [TJ we find that the effect 
of overcounting is an rs +0.7% shift in the measured value of Alr- This effect is small, 
however, and well within the theoretical error discussed above. 

^Certain corrections, such as those related the running of a, will induce effects that are suppressed relative 
to the one-loop result by a/71- In (M§ /Wg) ~ 5%; however, this is still a small shift, and a large class of these 
contributions can be resummed if better precision is required. 
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6 Conclusions 



We have presented a calculation of the complete O (a) radiative corrections to M0ller scat- 
tering, and have discussed their effect on the total cross section and polarization asymmetry 
measured in E-158, a low energy fixed target experiment being performed at SLAC. The 
virtual EW corrections have been previously computed in jllj : our recalculation serves as a 
check of this result, with which we find complete agreement. Our computation of the hard 
bremsstrahlung contributions with realistic experimental cuts is the first such result, and 
provides a complete theoretical prediction for comparison with the E-158 measurement. We 
find that hard bremsstrahlung effects induce a small ~ +4% shift in the measured value of 
Alr, the experimental cuts reduce the sensitivity of both the polarization asymmetry and 
total cross section to QED corrections. 

We have also reviewed the theoretical uncertainty in the prediction for Alr. The 
largest error arises from hadronic contributions to the 7Z vacuum polarization function W z , 
as discussed in : an improved dispersion relation analysis of low energy e + e~ — > hadrons 
data could reduce this value [22! ■ O ur results indicate that higher order QED corrections 
are likely to be small in comparison to the error on IF Z , and can be safely neglected in 
the E-158 analysis. The one-loop weak corrections are ~ 50%; although this effect seems 
dangerously large, it arises from the appearance of effects unsuppressed by the small electron 
vector coupling, a qualitatively new contribution which first enters at one loop. As no similar 
effect appears at higher orders, we conclude that the perturbative corrections to A LR are 
under control. 

The results obtained here are applicable to both low energy experiments such as E-158 
and future high energy colliders, and provide detailed Standard Model predictions to which 
measured values can be compared. We hope that they assist in maximizing the discovery 
potential of both current and future experimental programs. 
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Appendix 

We collect here the expressions for the various hard bremsstrahlung matrix elements \M. 
discussed in Eqs. l2~o1 and 1771 We first consider |.M^ a | 2 , which are computed with m e = 
everywhere; with the Mandelstam invariants presented in Eq. c^ (x) defined in Eq. and 
di defined in Eq. [T^l we can express them as 



\M h £L\ 2 = 8e 2 I (s 2 + s' 2 ) [tc 2 LL (t) d 2 d A + uc 2 LL (u) d 2 d 3 + t c 2 LL (tp) d x d 3 + u c 2 LL (up) d x d^ 
-2s [c LL (V) c LL (V) (s + t' + uj di + c LL (t) c LL (u) (s + t + u) d 2 } 
-2s [c LL (t'j c LL (u) (s + t + ujd 3 + c LL (t) c LL (u'j (V + t + u'^j d 4 } 
-2 [c LL (t) + c LL (u) } {c LL (t) + c LL (u)} s (t + u) (t + u) d x d 2 
-2 [c LL (V) + c LL (w)} [c LL (t) + c LL (u'j } s (V + uj (t + u'j d 3 d 4 
- {cll (u) c L l (t'j + c LL (Y) c LL (u'j + c LL (u) c LL (V ) } [tt (s + 8 



+ss (t — £ ) — ss (u + uj — uu ( s + s ) } d\d 3 

- [cll (t) c LL (u'j + c LL (uj c LL (t'j + cll (t) c LL (V)} [uu (s + s'j 

+ SS (u — U^j — SS (t + t'j — tt ( s + s ) } ^1^4 

- {cll (u) cll (*') + c LL (t) c LL (u) + c LL (t) c LL (Y ) } [uu (s + s') 
+ss (u — uj — ss (t + tj — tt ( s + s j I d 2 d 3 

- [cll (t) cll (uj + c LL (u) c LL (t) + c LL (u) c LL (V) } [tt (s + s 

+ SS (t — tj — SS (u + U^j — UU ( s + s ) } ^2^4 I ! 

\M&\ 2 = \M h I : a L (c LL (x)^c RR (x))\ 2 , 

\M h £ a R \ 2 = \M h jil\ 2 = 8e 2 \uc 2 LR (u)(t 2 + t 2 )d 2 d 3 + tc 2 LR (t)(u 2 + u 2 )d 2 d 4 



+t c \n (*') {s 2 + s ' 2 + t 2 ~ t 2 - 2t'u - 2t'u + 2st + 2ss + 2s t} d x d 3 
+u c 2 LR (u ^ [s 2 + s 2 + u 2 — u 2 — 2ut — 2ut + 2su + 2ss + 2s u\ d\d^ 
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+ [clr (£) c LR (V) {(u + u'^j (ss - tt + uu'^j + 2suu'} 
+clr (u) c lr (u'^j + ( ss ' - uu ' + tt"j + 2s££'}} did 2 
+ {clr (t) c LR (V) {(u + u'^j (ss - tt + uw'j + 2s'mt'} 
+clr (u) c lr (u'^j {{t + t') (ss - uu + ££') + 2s'tt'|| d 3 d 4 
+2c Lj r (u) c LR (u) t (s + u) (s + u) did 3 



+2c LR (t) c LR (V) u (s + t) (s + t) d x d A 
+2c LR (t) c LR (V) u(s + t) (s + t) d 2 d 3 

+2c LR (u) c LR (V) t(s + u') (s + uj d 2 d 4 j . (35) 

We now present |.M^' 6 | 2 , the mass terms which become important when — * (9 (ml); we 
find 

l-M^I 2 = 8e 2 ^ { - {c LL (Y) + c LL (u)} 2 s [s 2 + (t + uf] d\ + c\ R (t) u {t 

x (s — u — t) + su + 2ss +su}dl + c 2 LR (u^j t jii (s — u — t) 
+st + 2ss' + s't]d\ - {c LL (t) + c LL (u)} 2 s {s 2 + (t + uf] d\ 
+c 2 LR (£) u j£ (s — u — £ ) + su + 2ss + s -uj d\ + c 2 LR (u) t jii (s — u — t ) 
+st + 2ss' + s't] dj + 2 [c LL (u) + c LL (t) } 2 s 2 (u + t') d 2 . 
+2 {c LL (0 + c LL (u) } 2 s 2 (t + u) d\ | , 
l-M^I 2 = \M h / L (c LL (x)^c RR (x))\ 2 , 

\M h £ b R \ 2 = 8e 2 ^ | jc RR (Y) + Ci?R (m')} 2 s' (s + t + u'fdj - c\ R (t) u 

X |t + £ +«j+s?t+s«J^- c 2 LR (u^jt jii + t + It ) 
+s£ + sY} d 2 + {c LL (£) + c LL (u)} 2 s (s + t + uf d\ - c 2 LR (t) u 
x \t(s +t + uj+su + s «} c?2 — c 2 ^ (u) £ jit (s + £ + uj 
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+st + s't\ d\ + 2 c\ R (u) t (s + u)dj + 2 c 2 LR (V ) u (s + t') d\ 
+2 c\ R (t) u (s + t)dl + 2 <? LR (u) t (s + u) d\ j , 
= K^cuW-c^W)! 2 . (36) 

Finally, the modification of the total unpolarized cross section that results when t, u, t' ,u — > 
O {m 2 e ) is 



[s + u) (s + w') d 2 d 4 (s + t) (V + t') d 2 d 3 



unpoll "°e I ^2 ^2 

+ w'J fs' + Wj diC/3 + t'^j (s + t) d\d± 
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